Long Distance Entanglement Generation in 2D Networks 
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We consider 2D networks composed of nodes initially linked by two-qubit mixed states. In these 
networks we develop a global error correction scheme that can generate distance-independent en- 
tanglement from arbitrary network geometries using rank two states. By using this method and 
combining it with the concept of percolation we also show that the generation of long distance en- 
tanglement is possible with rank three states. Entanglement percolation and global error correction 
have different advantages depending on the given situation. To reveal the trade-off between them 
we consider their application on networks containing pure states. In doing so we find a range of 
pure-state schemes, each of which has applications in particular circumstances: For instance, we 
can identify a protocol for creating perfect entanglement between two distant nodes. However, this 
protocol can not generate a singlet between any two nodes. On the other hand, we can also construct 
schemes for creating entanglement between any nodes, but the corresponding entanglement fidelity 
is lower. 

PACS numbers: 03.67.-a, 03.67.Bg, 64.60.ah 



I. INTRODUCTION 

Quantum repeaters have shown great promise for cre- 
ating entangled states over long distances, which is es- 
sential for many quantum processing tasks [J-IH. The 
repeaters consist of a ID network of nodes linked by en- 
tangled states. By using local operations and classical 
communication (LOCC) an entangled state is generated 
between distant nodes, with the required number of ini- 
tial entangled states between each node scaling logarith- 
mically with the states separation distance @ . Once cre- 
ated the entangled state can be used for a range of tasks 
including quantum cryptography distributed quan- 

tum computing Q and quantum teleportation flfjj . The 
goal of these protocols is to generate a highly entangled 
state between two distant nodes and typically one uses 
the states fidelity, i.e. overlap with a singlet, as a measure 
of success in achieving this. 

However, the one dimensional nature of quantum re- 
peaters means that their application in higher dimen- 
sional networks will not take advantage of the higher 
connectivity By creating protocols for these higher di- 
mensional networks the required amount of initial entan- 
glement between the nodes can be reduced. This was first 
discussed for pure state networks and involved the use of 
'entanglement percolation' fldT[l6j . The procedure uses 
the 'procrustean method' to transform each pure state 
into a maximally entangled singlet with a finite proba- 
bility [13] ■ When this probability exceeds a geometry de- 
pendent threshold an infinite cluster is formed of linked 
nodes. Any two nodes in the cluster can then be trans- 
formed into a singlet using the process of entanglement 
swapping [l2|, EH • A singlet is generated between nodes 
with a probability that is independent of their separation. 
In particular, the required number of entangled states be- 
tween each node does not change with the separation, in 
contrast to quantum repeaters. 



Unfortunately, in reality every state will experience 
noise that causes it to lose its purity and become mixed. 
Therefore extending these ideas to mixed states is of vital 
importance for applications. In previous work we have 
extended entanglement percolation to networks that ini- 
tially consist of a particular type of mixed state [l9|, [2(| , 
that can occur when channels are subject to amplitude 
damping, and we have shown that these are the only two- 
qubit states that allow perfect singlets to be created. 

Other schemes have also been shown to generate highly 
entangled states from a regular 2D square network made 
of 'binary states' [!, [2l|. Here, if the amount of en- 
tanglement in the binary states exceeds a threshold, the 
network is transformed into an entangled state that can 
stretch over an arbitrary distance, while maintaining a 
non-maximal but constant entanglement. Our previous 
protocol [l^] , as well as the scheme described in Ref . [2l[ , 
requires states of rank two, or less, for long distance en- 
tanglement generation using constant resources. So far- 
no protocol has been able to transform a 2D network of 
full rank states into a highly entangled two qubit state 
with no dependence on the qubit separation, although 
this is possible in infinite 3D networks [22j]. 

In this paper we extend on these ideas to obtain a 
'global error correction' procedure that can be applied to 
quantum networks that have arbitrary geometry and are 
composed of binary states. By combining it with quan- 
tum state percolation we will show that entanglement 
can be efficiently generated over long distances by using a 
constant number of rank three states between each node. 
The final fidelity is independent of the distance. Pure 
states can also be used by the global error correction and 
entanglement percolation protocols. These can again be 
combined to create a range of altered schemes that allow 
for a continuous transition from a pure state global error 
correction method to entanglement percolation. We find 
that global error correction allows any two chosen nodes 
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FIG. I: (color online) a) A quantum network is shown com- 
posed of nodes (large circles) that contain qubits (black cir- 
cles). Each qubit is initially entangled to another (thin black 
line) and these two qubit states form the edges that are con- 
tained within each bond (thick grey lines), b) A 5x5 node 
square network with its dual network. The individual edges 
and qubits are not shown. However, the dual nodes are pre- 
sented as red, small circles connected by dashed lines, includ- 
ing the network exterior (dashed square with small circle), 
c) A 5x5 node triangular network with its honeycomb dual 
network. 



in the network to be linked by an entangled state, but 
this is at the expense of the state's fidelity. Entanglement 
percolation only allows two nodes from a known subset 
to be to be connected, but this is via a perfect singlet. 
If the criteria for entanglement percolation are not satis- 
fied we have to use a specific combined scheme to create 
states of maximal fidelity and the choice of this scheme 
is dependent on the initial states in the network. 

The paper is structured as follows. In Sec. HI we de- 
fine the network and develop a generalized form of global 
error correction. This scheme works on binary states 
and in Sec. [ill] we give a way of generating these binary 
states from rank three states. These capabilities are then 
combined to enable long distance entanglement distribu- 
tion with rank three states. Requirements for the scheme 
are calculated and compared with numerical results. In 
Sec. lIVI the schemes will be altered to work on pure states 
and this allows comparisons to be made between the pro- 
tocols. Finally we conclude in Sec. [V] 



II. GLOBAL ERROR CORRECTION 

A quantum network is constructed by sharing two- 
qubit entangled states between locations that we call 
'nodes'. All quantum operations must be local to these 



nodes but classical communication is allowed between 
them. The shared two-qubit states are referred to as 
'edges' and all of the edges between two nodes constitute 
a 'bond' (see Fig. [IJ. Ideally each edge will be a maxi- 
mally entangled singlet, however due to decoherence this 
situation is unrealistic and we must have a method to 
cope with errors. In this section we will devise a global er- 
ror correction method that performs local measurements 
and, by combining the results, information is obtained 
on the error locations. This information does not allow 
us to exactly locate the errors but, as we will see, does 
allow the errors to be subdued so entanglement can be 
created between distant nodes. Note that throughout we 
assume the measurements and operations act perfectly. 

We will consider the entire network as a mixture of 
pure states and analyze the protocol acting on each of 
these using the stabilizer formalism [23| . From these re- 
sults we can then determine the mixed state produced. 
To begin we construct the network's initial state us- 
ing the states of individual edges. Each of these initial 
two-qubit edge states can be transformed into a prob- 
abilistic mixture of Bell states [24|, via LOCC. This 
allows us to assume that each edge is in one of the 
four Bell states given by \ip ab ) = (X a Z b ® I)(|00) + 
1 1 1) ) / , where a, b € {0, 1}. A stabilizer generated by 
{{-l) a Z <g> Z, {-l) b X <g> X) can then describe each Bell 
state and we label the probability of that state as p a b- 
For a network that contains Ne of these identical edges, 
and Nq = 2Ne qubits, the whole state will have the form 

N E 

Pi = YlP^biH'a^) (^a,b,\ 

{6|}{Oi}i=X 

where each edge, i, has parameters £ {0, 1} and 6j G 
{0,1}. The summations are over all values that these can 
take and give a mixture of pure states 



(2) 



that occur with probability 



P{{a i },({b i }) = l 



(3) 



The error model we consider has independent bit-flip and 
phase-flip errors. For this case the pure states that con- 
tribute to the mixture can be thought of as networks 
of error free singlet states, IV'oo); on which errors ran- 
domly occur. A bit-flip error that causes a — > 1 oc- 
curs with probability p x . Similarly, phase- flip errors 
cause b — > 1 and have a probability of p z . This sit- 
uation is equivalent to having edges with the values 

POO = 0--Px)0-~Pz), P01 = {l-Px)Pz, PlQ = (l-p z )p x , 
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and pn = PxPz and it simplifies the probability of each 
state P{ ai },{bi} to become 

P({a. ( } , ({b,}) = p^ x (1 - Px ) Ne - Nx x 

p»Hl- Pz ) N *- N *, (4) 

with Nx and Nz giving the number of bit- flip and phase- 
flip errors on the state, respectively, 

i=l 
N E 

N z =Y J h- (5) 

i=l 

Now we apply the global error correction procedure to 
one of these pure states. Each of which can be described 
using a stabilizer that is formed by the union of Bell 
state stabilizers. The procedure extracts information by 
performing Z <£> Z measurements locally on nodes along 
a closed path of edges, as illustrated in Fig. Such 
paths arc referred to as 'loops' and contain edges from 
the set loop = {h,l2, ■■■■In l }, where the give the 

index of each of the Nl edges forming the loop. These 
measurements take —1 and 1 as possible values. The 
product of the measurement results around a loop is la- 
belled (— l) L i°°p. Here Li oop gives the bit-flip error par- 
ity around the loop, L loop = (jL, ie i oop a i) m °d2. Each 
value of Li 00 p gives us information on the errors that we 
can then use. As the measurements are made the sta- 
bilizer follows the rules given in Ref. [2^]. After each 
measurement the stabilizer generator is manipulated so 
that at most one element anti-commutes with the mea- 
surement operator and the measurement operator, mul- 
tiplied by the outcome, is then substituted in place of 
the anti-commuting element. If the operator does not 
anti-commute with any elements in the stabilizer then 
nothing is changed. 

We obtain the measurement outcomes, (— l) mj (rrij £ 
{0, 1}), for the operator Z <£> Z between every pairing of 
qubits in each node. By making these measurements we 

can calculate Li oop = (j2i<zi oop m ij rnod2 for all of the 
smallest possible loops in the network, which we refer to 
as plaquettes, and can then deduce Li oop for any loop. 
We do not need to perform the actual measurement be- 
tween every qubit pair in a node because some of the 
Z ® Z operators will be products of previously measured 
Z <g> Z operators. From the state P{ ai },{bi} these mea- 
surements result in a state given by a stabilizer that has 
generators (-l) Nz (X® N Q), (-l) a '(Z®Z) for every edge, 
i, and operators (— l) mj (Z ® Z) for every measurement 
3- 

All of the plaquettes in the network can then be as- 
signed their values of Li oop - Values of one can be consid- 
ered as defects on nodes in the dual network (see Fig. [2]) . 




FIG. 2: (color online) One possible pure state network 
P{a },{b } from the ensemble given by Eq.[T] Each node (white 
circles) contains qubits (black dots) that are part of two qubit 
entangled states (thin lines) and edges are located in bonds 
(thick grey lines). For a = the edge line is solid and if 
a = 1 it is dashed. Z®Z measurement operations are carried 
out within each node (black ovals) around one plaquette. By 
doing this a value of Li oop can be associated with each plaque- 
tte. The dual network (red dashed lines) is shown with these 
values contained within each dual node (red circle). The dual 
vertex representing the exterior is represented by a dashed 
box. 



These defects occur in pairs. A chain of errors, <2j = 1, 
on the quantum network's edges separate the dual node 
defects in each pair (see Fig|3k). These error affected 
edges must cross a path of dual network edges that con- 
nect the two defects. Hence, we can obtain information 
on the location of bit-flip errors by pairing the defects, 
in a matching problem. We assume that the length of 
the paths is minimal, since longer paths require more er- 
rors and this is unlikely when p x is small. This gives us 
a minimum weight matching problem to be solved pl| 
in an analogous way to the correction methods used in 
surface codes |26h28| . By doing so we have found the 
most likely defect pairs to be linked by a path of errors. 
The actual path is assumed to be the one of shortest 
distance, i.e. smallest number of edges, linking the de- 
fects on the dual network (see Figf3^). Sometimes more 
than one such path exists and we would wish to mini- 
mize the number of nodes between the true path and the 
assumed path. This can be easily done on regular net- 
works, by using paths that approximately follow straight 
lines, but for simplicity in the general case it is sufficient 
to choose one of the shortest paths at random. One qubit 
in each network edge that lies along the assumed paths 
is then acted on with a X gate. We can actually avoid 
performing these X gates until the very end when we may 
have fewer qubits to apply gates on. Here we include the 
gates at this stage for clarity. When these gates are ap- 
plied to a qubit we also change the value of rrij to 1 — rrij 
for each measurement pairing that includes those qubits 
(see FigODa). This gives rrij = 1 for an even number of 
times around every loop and by acting around every loop 
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FIG. 3: (color online) a) A single edged network is shown 
where three edge errors have occurred (dashed lines). The lo- 
cal Z ® Z measurement outcomes are displayed inside boxes, 
with '±' referring to a ( — l) mj = ±1 eigenstate. Around each 
plaquette the values of Li oop axe then calculated and shown 
on the dual network nodes (red outlined circles linked by dot- 
ted lines). Values of 1 are defects and are matched together 
along a path (blue arrows), b) Edge errors are assumed to 
lie across the matching path and X gates are applied to one 
qubit within each of these edges (black arrows). This 'flips' 
(+ -H> — ) some of the local measurement outcomes. Each 
qubit can then be put into two groups (labelled '0' and '1') 
using the measurement outcomes ('-' means the groups need 
to be different for the two qubits) and qubits that are in a 
network edge get put in the same group. To assign the val- 
ues we pick an initial qubit and value at random and then 
assign values in all remaining qubits. c) X gates are applied 
to one group (nodes containing white X's) and for each pure 
ensemble state a GHZ state is created with patches (shaded 
red area) of bit-flip errors. These occur where the actual edge 
error path differed from the assumed path. 



with X gates we can force rrij = for every local qubit 
pairing. A procedure to accomplish this task is given 
in Figj3jb-c). Once this is done the resulting stabilizer 
will have generators (-1)^^®^, (~\)<Z <g> Z for ev- 
ery edge, i, and Z '(g)Z from every local pairing of qubits. 
The a, values have been changed to a[ after performing 
the X gates. These a! i values are unknown, however all 
of the dual edges that cross edges with a\ = \ must form 
loops on the dual network. This stabilizer describes a 
pure GHZ state that has had Nz Z operators acting on 
it and X operators acting on all the qubits within some 
nodes. These affected nodes are within loops on the dual 
network and appear as patches (see Fig. [4|. The X gates 
that remain when the procedure is run on an initial con- 
figuration {a,i} arc described by an operator Xp({a,i}), 
that applies X gates on the qubits that lie within patches. 
The whole procedure will result in a mixture of these pure 
states, 

Xp({a l })Z Nz PGHZ Z Nz X P ({a l }), (6) 

where Z acts on an arbitrary qubit and P({at} | {rrij}) 
is the probability of the initial error configuration being 
{cii} given the results of the measurements. 

Larger patches require more errors and are less likely 
to occur. Hence, for low values of p x the most prob- 
able case is that finite sized patches, containing a small 
number of nodes, occur randomly. We quantify the prob- 
ability of a random node existing in a patch as Px , which 
has the same value throughout the network. For larger 
values of p x the likely number of errors increases so that 
it becomes impossible to match the defects up correctly 
and accurately guess the paths. Within infinite networks 
there exists a critical value of p x below which patches al- 
most certainly exist and long distance entanglement can 
be generated. Above the critical value this is not possible 
and Px = 1/2 for the most likely pure states. 

When we are highly likely to form patches the resulting 
mixed state can be cut back to an impure GHZ state, 
with any number of qubits up to Nq, by removing qubits 
with X measurements and, depending on their outcomes, 
applying a correcting Z operation on one of the qubits to 
be kept. The final fidelity, i.e. overlap with a pure GHZ 
state, is dependent on p x and decreases for more qubits 
in the resulting state. It is however still independent 
of the distance between the qubits. In particular, we 
can remove all of the qubits except two leaving a two- 
qubit entangled state spread between two nodes. For 
each pure state contributing to the mixture in Eq. Q 
after removing the qubits we obtain an ideal singlet when 
the remaining qubits are either both in an error affected 
patch or outside of them (see Fig. The pure state 
P{a,i},{bi} is transformed to the state p c {{bi} , that 
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FIG. 4: (color online) Regular square network formed by sin- 
gle edged bonds. The dual lattice is shown in red/dashed 
and includes the error syndromes on the dual vertices (red, 
small circles; Li oop = l-filled, Li oop — 0-outlined). The edge 
errors causing these plaquette errors are shown by thick blue 
lines. The error syndromes are matched together (black ar- 
rows) and this can be used to give the most likely edge error 
locations for the recorded measurement syndrome. The grey 
areas between the assumed error locations and actual edge er- 
rors identifies groups of nodes. For this particular pure state 
in the ensemble we will obtain a GHZ state with X errors 
present on the qubits within these nodes. In practise we do 
not know what the initial state was and hence do not know 
the location and size of the grey patches, but for the regions 
to become larger more errors are required and hence are less 
likely. 



has the form 



Pc ={P 2 X + (1 - Px?) 



l^ooX^ool 



2Px(l Px) {1 + (1 2Pz)JVE W )frM 



(P| + (1-P x ) 



(l-(l-2p,)"*) 



\^oi}{ipoi\ 



(1- (I - 2r> ) Ne ) 
+ 2P X (1-P X ) [ - 1 2 Pz) V i)frM, (7) 

and the final mixed state is given by the weighted sum 
of these, 

pf= E ^(i-f/ B ^%}iw)x 

Pc {{bi} ,{oi}). (8) 

Since Px is independent of the state separation we see 
that as long as there are no phase errors present, i.e. 
p z = 0, the state produced has a fidelity that is inde- 
pendent of the node separation. These states that do 




FIG. 5: (color online) Here a 20x20 network has had its bonds 
removed with a probability of 0.05 and then X errors intro- 
duced with a probability of 0.05. We apply our global er- 
ror correction procedure and eventually this generates a GHZ 
state with patches of nodes that exhibit unwanted, extra bit- 
flips (red, filled circles). If we try to link a node from within 
a patch to outside of them (dashed line) the final state two- 
qubit state exhibits a bit-flip error. Examples are also given 
of qubit pairings that would not exhibit errors (bold black 
lines). 



not experience phase errors are called binary states 
Our global error correction procedure is a generalization 
of the one given in pll ] but does not involve teleporting 
qubits and can be easily applied to any network geometry. 
Later we will see that we can randomly transform states 
into binary states, obtaining them on random bonds but 
also destroying all of the entanglement in the remaining 
bonds. Our scheme can still be applied in these situa- 
tions and in Fig. [5] it has been applied on a contribut- 
ing term in the largest cluster of a square network that 
was missing such bonds. We will discuss this in more 
detail in the next section (see Sec. IIII| . Furthermore, 
in regular networks with all bonds intact we performed 
simulations to reveal the critical threshold of p x below 
which long distance entanglement distribution is always 
possible. This was done for both square and triangular 
networks (see Fig. QJb-c)). In both cases we see the final 
states fidelity decrease suddenly after critical p x values 
(see Fig. . These transitions become more prominent 
for larger networks and later we will attach values to the 
thresholds. 

If we have both bit-flip and phase errors present the 
initial edges will be rank four Werner states [29[ . To ap- 
ply the procedure in this case we can use error correction 
codes to subdue the phase- flip errors 0, [3(| [3l|. All 
of the qubits are then replaced with encoded qubits and 
all operations are replaced with their encoded versions. 
This technique was looked at in pll ] for a square network 
using a majority voting CSS code. In this case for bonds 
containing 2t + 1 edges the phase error probability on 
each encoded edge can be suppressed to w 
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FIG. 6: The fidelities are plotted from simulations of global 
error correction when run on complete square and triangu- 
lar networks to produce a two qubit entangled state between 
two randomly chosen qubits. In the square case (solid lines) 
we see a lower threshold compared to the triangular network 
(dashed lines). We also included the approximate threshold 
values (dotted vertical lines) that are calculated in Sec. Mil 
We used networks of 10x10 nodes and 20x20 nodes, with the 
larger networks presented by thicker lines, and can see that 
the threshold becomes more pronounced when the network 
gets larger. 



To maintain a constant phase error probability on the 
final state the number of edges in each bond only needs 
to increase logarithmically with the size of the network. 
Unfortunately, this operation also increases the chance 
of a bit flip on the encoded edge to sa (2t + l)p x , which 
means there is a point where the chance of a bit-flip ex- 
ceeds the network's critical value, patches are no longer 
small and the procedure fails. When p x is small this still 
enables the p z error probability to be substantially re- 
duced. For rank two, binary states, the number of edges 
required in each bond does not need to scale with dis- 
tance and yet for rank four states there is a logarithmic 
scaling. In the following we are interested in extending 
the constant scaling with states that are of rank three. 



III. ENTANGLEMENT DISTRIBUTION WITH 
RANK THREE STATES 

Using the global error correction procedure we can 
cause the edge errors in every pure state in an ensem- 
ble to condense into patches on a large GHZ state in 
such a way that a highly entangled two-qubit state can 
be generated between two distant nodes. This requires 
the initial network nodes to be linked by sufficiently en- 
tangled binary states. However, in the following we will 
show that such states can be probabilistically generated 
from a larger class of rank three states using LOCC. This 
widens the class of states in a 2D network that allow for 
long distance entanglement distribution. 

It is impossible for a rank two, binary state to be cre- 



2D Square 0.5 

2D Triangular 2sin(7r/18) w 0.347 

2D Honeycomb 1 - 2sin(7r/18) g 0.653 

TABLE I: Threshold probabilities for various regular network 
geometries [53. l33|. 



ated using anything other than states that arc rank three 
or less. We give a short proof of this fact in appendix |A"1 
If we have two rank three states that can be probabilis- 
tically transformed by LOCC into states of the form 

p(X,u) = A|VooX^oo| + *#oi)(V'oi| + (1 — A — t/)|01><01|, 

(9) 

with < A, v < 1 and A + v < 1, then they can be 
probabilistically distilled into a binary state 

Pb (p' x ) = (i -pDhMftM + {p' x )Ww){<t>w\, (io) 



where 



P., 



1 



2{\ + v y 



(ii) 



The states given by Eq. © can actually form (up to 
local unitaries) when both qubits of a maximally entan- 
gled state, IV'io); pass through phase flip and amplitude 
damping channels. Having two of these states the trans- 
formation to a binary state is done by performing a pro- 
tocol discussed in Ref. pjl, [2(|, where it was called a 
'pure state conversion measurement' (PCM), followed by 
local Hadamard operations. The PCM involves perform- 
ing two C-NOT gates locally, with one entangled state's 
qubits acting as the target qubits. These target qubits 
are then measured in the computational basis. If we find 
both qubits to be in the state |1) it has succeeded, which 
happens with a probability of (A + v) 2 /2. If it fails then 
we destroy the entanglement and lose both edges used. 
A graph of p' x and the probability of succeeding is shown 
in Fig. [7] From this it can be clearly seen that to min- 
imize the resulting p' x we would like |A — v\ to be max- 
imal. The maximal success probability with two states 
is 1/2. If we have m initial states in a bond we can 
generate a binary state with finite probability exceeding 
P c = 1-(1- ((A + ^) 2 /2)) Lm/2j , which can be arbitrarily 
close to one given enough edges. Again, if this fails then 
all of the edges are lost. 

By using this distillation procedure it is possible to 
transform an initial network composed of rank three 
states into a network of binary states, where the bonds 
are missing with a certain probability. This is the case 
for the network shown in Fig. [5] All of the nodes that 
are connected by a path of entangled states are then said 
to be in a cluster. For entanglement to be generated be- 
tween two nodes they must both be in the same cluster 
and this links with results from bond percolation the- 
ory (32| . For a square network of infinite extent, if the 
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FIG. 7: a) p' x coefficient for the binary state created using 
a PCM. b) The probability, P c , of generating a binary state 
from two initial states, i.e. m = 2, of the type given by Eq.[9] 



success probability of transforming each bond into a bi- 
nary state is larger than 1/2 then a cluster that contains 
infinite nodes is almost certain to exist. This is the 'per- 
colation threshold' and values for further geometries are 
given in Tab. [I] The probability of a node being inside 
the infinite cluster is given by the 'percolation proba- 
bility', <j>(P c )- The values of these together with ip(P c ), 
the probability that a bond exists and is linking nodes 
in the infinite cluster . can be calculated numerically for 
different networks 34 , HI | • For large finite networks this 
threshold phenomenon still exists but the transition be- 
comes smoother. In this case, when the critical value 
is exceeded we obtain a large cluster that contains a 
majority of the nodes. On this cluster we can then ap- 
ply our general global correction procedure to obtain a 
highly entangled state between two nodes. This requires 
that the binary states present in the largest cluster are 
sufficiently entangled for global error correction to suc- 
ceed. We use entropic arguments to give an approxi- 
mate threshold for a finite square network, with sides 
consisting of L nodes. The cluster will contain an aver- 
age of (Ne) = 2L(L — l)ij>(P c ) binary states, with each of 



these introducing a Shannon entropy H(p' x ). The mea- 
surements then extract one bit of information from every 
plaquette in the cluster. The number of these plaquettes 
on average is given by (Np) = 2 + (Ne) — (Nn), with 
(Nn) = L 2 <j)(px) being the average number of nodes in 
the cluster. We require (Np) > (Ne) H(p' x ) for enough 
information to be gathered and this gives a bound for 
p' x when we solve H(p' x ) = (Np) / (Ne)- These bounds 
define a region for P c and p' x within which long distance 
entanglement distribution can be achieved. In the infi- 
nite network case, N — > oo, we have ip(P c ) = Pc4>(Pc) 
and can calculate that 



(N P ) I (N E ) = 1 + 



N- 



2 - L 2 4>(P C ) 
2L(L - l)P c cj>(P c ) 

1 - {2P C )~ X . 



(12) 



This defines the region by the relation (1 — H(p' x ))P c > 
1/2 and the boundary is shown in Fig. [5J Our procedure 
is also related to methods developed to cope with qubit 
loss in surface codes (3(| and similar critical regions have 
been found in these situations. For the case where P c = 1 
the critical value of p' x in an infinite square network is 
then given by H(p' x ) = 1/2. This can be solved to yield a 
threshold of p' x = 0.11. Further analysis for the complete 
square case is performed in Ref.[2l[ and by relating the 
scheme to surface code error correction [281 ] a threshold 
of p' x = 0.1094 is found. Similarly, in a complete, infi- 
nite triangular network we have H(p' x ) = (Np) / (Ne) = 
l + (2-L 2 )/( J L-l)(3L-l) -> 2/3 yielding a threshold of 
p' x = 0.17. Both of these thresholds are shown on Fig. [6] 
In this figure it can be seen that the higher connectivity 
present in the triangular network does provide for higher 
fidelity and a lower threshold than a square network. The 
simulations were run for finite networks so a perfect tran- 
sition at the critical threshold does not occur. However, 
as the networks become larger the fidelity will approach 
1/2 for our critical values. Note that these critical values 
are calculated for the ideal method and will be slightly 
larger than the simulated threshold due to the random 
choice of paths. 

We have performed numerical simulations of the 
scheme, when edges are removed randomly, for a 25x25 
square network. The final singlet state fidelity for differ- 
ent distillation parameters was calculated and the results 
are shown in Fig. [8] Our boundary between the values al- 
lowing long distance entanglement distribution and those 
that do not is also shown. This gives a good separation 
between the parameters producing singlets and those re- 
sulting in separable states. There is no sharp transition 
between the low and high fidelity region which is due to 
the finite network size. The limiting cases of pure state 
percolation and binary state global error correction are 
given by P c > 0.5, p x = and P c = 1, p' x < 0.11, respec- 
tively. 
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FIG. 8: (color online) a) The probability of a node existing in 
the largest cluster, 4>. b) Based on a 25x25 square network, the 
fidelity F between two randomly chosen nodes in the largest 
cluster of connected nodes is shown depending on the bond 
conversion probability P c and X-error probability p' x . The 
boundary for values that allow long distance entanglement 
distribution in infinite networks is shown as a thick black line. 
Pure states can be probabilistically transformed into binary 
states with parameters along the dashed lines. These lines 
are for the pure states a = 0.75, 0.76 • ■ • 0.85 from lowest to 
highest p' x values, respectively. 



IV. TRADEOFFS BETWEEN DIFFERENT 
STRATEGIES ON PURE STATE NETWORKS 

It is interesting to look at the behavior of entangle- 
ment percolation and global error correction when cre- 
ating long distance entanglement from the same initial 
network. All of these procedures require a network of 
identical initial states. If the bonds in the network con- 
tain enough entanglement we produce a highly entangled 
state between two nodes that are a long distance apart. 
The two nodes that are connected are out of a random set 
of nodes. Different methods may have benefits such as re- 
quiring less entanglement in the initial bonds, generating 
higher entangled states between two nodes or allowing a 
greater proportion of nodes to be linked by the entangled 
state. Pure states are the simplest case where this com- 
parison can be considered so we will look at this case. By 
randomly generating states between nodes from the ini- 
tial states and then using global error correction we have 
found a method that works on networks containing rank 
three states. By changing the states generated we can 
apply this scheme on networks of pure states, together 
with global error correction and entanglement percola- 



tion. We will see that there is a wide choice of states 
to generate and each of these may provide further bene- 
fits compared to entanglement percolation or global error 
correction. 

For pure non-maximally entangled states entanglement 
percolation is known to succeed in creating a long dis- 
tance singlet state deterministically [111 ]. This requires a 
geometry dependent amount of entanglement, related to 
the percolation thresholds (see Tab. HJ. The pure states 
are characterized by their Schmidt decomposition 



| a) = Va\00) + VI 



(13) 



with a being the largest Schmidt coefficient, 1 > a > 1/2. 
Any pure state can be put into this form via local unitary 
operations and a can be considered as a measure of the 
state's entanglement. 

Apart from basic entanglement percolation other tech- 
niques have also been proposed in pure state networks 
with smaller initial entanglement requirements [llTtl3l 
[l5j | . These include the global error correction on a square 
network with no missing edges, as discussed in pl| . and 
it was shown that a lower amount of entanglement was 
required to succeed at the expense of the final states' 
fidelity. By randomly generating different states before 
global error correction it becomes possible to explore the 
transition between global error correction and entangle- 
ment percolation. Here we will again choose to create 
binary states by probabilistically converting pure states 
into the binary states. When we start with pure states 
there are a number of ways to perform this operation. 
These create the binary states with higher probability of 
success or more entanglement. To perform the conver- 
sion each pure state is transformed, using LOCC, into 
another pure state 



\a) = Ve7|00) + Vl - a'|ll), 



(14) 



that has higher entanglement with the optimal probabil- 
ity of succeeding P c = (1 — a)/(l — a'). This operation 
involves a local measurement at one qubit. The result 
is then communicated to the other qubit where a uni- 
tary is performed. This transformation operation and 
probability come from 'Majorization' results [37[. These 
states can then be 'twirled' @, HH into binary states, 
^) 2 /2. 



withj4 = (vff - VT — a') 2 /2. These actions probabilis- 
tically generate binary states at random bonds, allowing 
global error correction to be applied. The value of P c 
can be adjusted between and 1. If P c is larger than 
the percolation threshold then an infinite cluster, of bi- 
nary states pb(p' x ), is formed and a node is a member 
of this cluster with probability <j)(P c ). Any two nodes in 
the cluster can then be linked by global error correction. 
Two nodes are in the cluster with a probability of (j>(P c ) 2 
and they can be turned into a binary state that has fi- 
delity F after global error correction. In Fig. [8] the paths 
of (1 — P c , p' x ) are shown for different values of a as a 
dashed line. Each point on these lines correspond to a 
binary state that can be obtained from the initial pure 
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FIG. 9: (color online) Starting from pure states we can create 
a range of binary states (dashed lines in Fig. [8}. For each of 
these the largest cluster contains a different proportion (f> of 
total nodes and any two of these can be linked using global 
error correction to leave a two qubit entangled state that has 
fidelity F. The different values of F and <fi are given for each 
initial pure state. Shown are plots corresponding to eleven 
initial pure states with a = 0.75, 0.76, 0.77, . . . , 0.85 (top to 
bottom). 



state. The relationship between the available fidelities, 
F, and fraction of nodes in the cluster is given in Fig. |H1 
Each line relates to a different initial pure state. For each 
dashed line representing an initial state in Fig.|5]the only 
points allowing long distance entanglement in the infinite 
case lie below the critical boundary, shown in Fig. [3] 

If we are required to generate long distance entangled 
states with the highest possible amount of entanglement 
then, given that the threshold for entanglement perco- 
lation is satisfied, the preferred method is to randomly 
obtain perfect singlets and then use entanglement swap- 
ping. This generates a perfect singlet, however in doing 
so we reduce the number of nodes that are available to be 
linked. When the threshold is not satisfied we must prob- 
abilistically transform the initial states to binary states 
with the largest 1 — Pq that still lie within the critical re- 
gion. By doing this one highly entangled state is created 
over a long distance. Note that the initial state of the 
edges needs to be known in order to decide on the state 
to percolate. We are unable to choose the nodes before 
running the protocol, but we still know which nodes are 
linked at the end. If we are required to link two partic- 
ular nodes, that have been chosen before the protocol is 
run, we must maximize 4> and this is achieved by using 
global error correction. The choice between maximizing 
the fidelity or (f> is dependent on the intended use for 
the entangled state. If one entangled state between any 
two nodes in a set is sufficient it may be beneficial to 
maximize <p. In order to judge the schemes a figure of 
merit can be constructed from (j) and the fidelity. This 
can then be maximized to reveal the preferred scheme. 
For example a suitable figure of merit in the case where 
two specific nodes need to be entangled would be F(j) 2 
and this is maximized for global error correction. 



V. CONCLUSION 

Previously, entangled states could only be generated 
over an arbitrarily long distance from rank two states, or 
less, in regular 2D networks [TTJ, [l9l - |2lj ]. Here we have 
devised a global error correction scheme that enables 
a highly entangled state to be generated over an arbi- 
trary distance. This procedure can be applied to any 
network composed of binary states and allows the cre- 
ation of multi-qubit GHZ states. By combining this with 
entanglement percolation we have extended this ability 
to any network and a class of rank three states. These 
states would result from a combined amplitude damping 
and phase error channel. For networks composed of these 
rank three states it becomes possible to produce highly 
entangled states over arbitrary distances, with constant 
resources between the nodes. Although this is still a re- 
stricted case it is one step closer towards 2D networks 
composed of full rank states and allows us to deal with 
two important kinds of noise. Initially there must be 
a sufficient amount of entanglement between the nodes 
for the scheme to succeed and we have provided an en- 
tropic estimate for the requirements in percolated net- 
works, based on the estimate in Ref. 2l|. We have 
numerically calculated the resulting fidelities and the re- 
sults agree with the threshold estimate. The behavior 
of both entanglement percolation and global error cor- 
rection was previously investigated separately but here 
we have provided a method to connect these possibilities 
and in doing so revealed the preferred methods for dif- 
ferent priorities. It is still an interesting question to ask 
whether long distance entanglement distribution can ex- 
ist in a 2D network of Werner states even when the gates 
are assumed to be perfect. 

This research was supported by the ESF via Euro- 
QUAM (EPSRC grant EP/E041612/1). 



Appendix A: Impossibility to generate binary states 
from rank four states 



Binary states can not be generated from two-qubit full 
rank states and we prove this by contradiction. If there 
were to exist some local procedure that produced some 
binary states from a full rank state then, by acting af- 
ter with local projective measurements on all qubits ex- 
cept those in one binary state, we can easily construct a 
method to generate a single binary state and pure sepa- 
rable state. This local procedure can be described by an 
operator Ma <£) Mb, where Ma and Mb act locally. The 
dimension of the entire Hilbert space is D. They must 
also satisfy M a m\ and N B N B < 1. The initial mixed 
state is 



D 
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with pi > 0, J2iPi = 1 an( i I*) £ ^4 ® is a complete 
basis. The procedure can then convert this initial state 
into 

D 

^2p t M A (g) N B \i)(i\M\ ®Nl= pp bm ® Psep- 

i=l 

p se p is a pure separable state, pun is the two-qubit binary 
state and p a non-zero probability. Since the initial state 



was full rank this property requires Ma <E> Mb to be rank 
two. This is only the case when either Ma or Mb have 
rank one. Any operator, Ma ® -Mb, with this property 
would not generate an entangled state. Hence, a binary 
state can not be generated using LOCC from a full rank 
state and this includes the case of finite full rank two- 
qubit states. 
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